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Abstract 

We revisit the issues on the thermodynamic property of stellar self-gravitating 
system arising from Tsallis' non-extensive entropy. Previous papers (Physica A 
307 (2002) 185; ibid. (2003) in press (cond-mat/0204315)) have revealed that the 
extremum-state of Tsallis entropy so-called stellar polytrope, has consistent thermo- 
dynamic structure, which predicts the thermodynamic instability due to the nega- 
tive specific heat. However, their analyses heavily relies on the old Tsallis formalism 
using standard linear mean values. In this paper, extending our previous study, we 
focus on the quasi-equilibrium structure based on the standard framework by means 
of the normalized ^-expectation values. It then turns out that the new extremum- 
state of Tsallis entropy essentially remains unchanged from the previous result, i.e., 
the stellar quasi-equilibrium distribution can be described by the stellar polytrope. 
While the thermodynamic stability for a system confined in an adiabatic wall com- 
pletely agrees with the previous study and thereby the stability /instability criterion 
remains unchanged, the stability analysis reveals a new equilibrium property for the 
system surrounded by a thermal bath. In any case, the stability /instability criteria 
are consistently explained from the presence of negative specific heat and within 
the formalism, the stellar polytrope is characterized as a plausible non-extensive 
meta-equilibrium state. 
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1 Introduction 



In many astrophysical problems involving self-gravitating many-body system, 
thermodynamic or statistical mechanical treatment usually loses its power 
due to the peculiar feature of long-range attractive force. In some restricted 
cases with the long-term evolution of stellar self-gravitating systems, however, 
thermodynamic stability criterion or statistical mechanical analysis recovers 
the physical relevance and plays an essential role in predicting the fate of 
such system. In fact, it is well-known that late-time stellar dynamical evolu- 
tion of the globular clusters as a real astrophysical system is driven by so- 
called gravothermal catastrophe, i.e., thermodynamic instability arising from 
the negative specific heat, which is widely accepted as a fundamental astro- 
physical process [1,2,3]. Historically, the gravothermal catastrophe has been 
investigated in detail considering the very idealized situation, i.e., stellar self- 
gravitating system consisting of many particles confined within a cavity of 
hard sphere [4,5,6,7,8]. In particular, special attention to the statistical me- 
chanical approach using the Boltzmann-Gibbs entropy has been paid. 

Recently, this classic issue was re-analyzed by present authors on the basis 
of the non-extensive thermostatistics with Tsallis' generalized entropy [9] (for 
comprehensive review of Tsallis' non-extensive formalism and its application 
to the other subject of physics, see Refs. [10,11]). In contrast to the Boltzmann- 
Gibbs entropy, the quasi-equilibrium distribution 1 

characterized by the Tsallis entropy can be reduced to the stellar polytropic 
system [12, 13]. Evaluating the second variation of entropy around the quasi- 
equilibrium state, we have developed the stability analysis and discussed the 
criterion for onset of gravothermal catastrophe (Ref.[14], hereafter paper I). 
Further, to clarify the origin of this instability, thermodynamic properties has 
been investigated in detail calculating the specific heat of the stellar polytrope 
(Ref. [15], hereafter paper II). The most noticeable thing in their papers is that 
the existence of thermodynamic instability indicated from the second variation 
of entropy or free-energy is completely explained from the presence of negative 
specific heat. As a consequence, the gravothermal instability appears at the 
polytrope indices n > 5 for a system confined in an adiabatic wall and at 
n > 3 for a system surrounded by a thermal bath. 

While the results in papers I and II are indeed satisfactory and the phys- 
ical interpretation of the instability is fully consistent with previous early 

sakagamiOphys.h. kyoto-u.ac.jp (Masa-aki Sakagami). 

1 Here and in what follows, we use the term quasi-equilibrium rather than the (ther- 
mal) equilibrium, since there is no strict thermal equilibrium in a self-gravitating 
system. The quasi-equilibrium state means that the system is, at least, stable in a 
dynamical equilibrium state. 
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works, from a standard viewpoint of the Tsallis non-extensive formalism, sev- 
eral important issues still remain unresolved. Among these, the most crucial 
problem is the choice of the statistical average in non-extensive thermostatis- 
tics. The previous analyses have been all investigated utilizing the old Tsallis 
formalism with the standard linear mean values, however, a more sophisti- 
cated framework by means of the normalized g-values has been recently pre- 
sented [16,17]. As several authors advocated, the analysis using normalized 
g-values is thought to be essential, since the undesirable divergences in some 
physical systems can be eliminated safely when introducing the normalized 
g-values. Further, non-uniqueness of the Boltzmann-Gibbs theory has been 
shown using the normalized q- values [18]. Of course, this does not imply that 
all the analyses with standard linear means or un-normalized q- values lose the 
physical significance 2 , however, in order to pursue the physical reality of the 
non-extensive thermodynamics, thermodynamic structure of the stellar self- 
gravitating system still needs to be investigated and the relation between the 
old and the new formalism must be clarified. 

This paper especially focuses on this matter starting from the construction 
of the extremum-entropy state by means of the normalized g-values. Within 
a mean- field treatment, we investigate the thermodynamic property of the 
quasi-equilibrium distribution surrounded by an adiabatic and a thermally 
conducting wall. In particular, we discuss the existence or the absence of ther- 
modynamic instability evaluating the specific heat of the quasi-equilibrium 
systems. Further, in order to check the consistency of the thermodynamic 
structure, the stability/instability criteria obtained from the thermodynamic 
property are re-analyzed from the second variation of entropy and free-energy. 
We found that the extremum state of the Tsallis entropy essentially remains 
unchanged and is described by the stellar polytrope. While the stability of the 
system in an adiabatic wall exactly coincides with the results in paper I, the 
new formalism using the normalized g-values reveals a new thermodynamic 
structure for a stellar system surrounded by a thermal bath. This point will 
be discussed in detail comparing it with the previous results. 

The paper is organized as follows. In Section 2, employing the standard Tsallis 
formalism, we revisit the issue on the most probable state of stellar quasi- 
equilibrium distribution on the basis of maximum entropy principle. While 
the resultant extremum state of Tsallis entropy reduces to the same stellar 
polytropic distribution as previously found, the new equilibrium distribution 
has several distinct structures. Taking fully account of this fact, in Section 
3, the thermodynamic properties of stellar polytropic system are investigated 
in detail. Thermodynamic temperature in stellar system is identified through 
the modified Clausius relation, which is indeed consistent with the recent claim 



2 Indeed, even the old formalism consistently recovers the standard Legendre trans- 
form structure leading to the usual thermodynamic relations [19,20]. 
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based on the thermodynamic zero-th law. Then, we evaluate the specific heat 
and discuss the existence or absence of thermodynamic instability. In Section 
4, the stability/instability criteria are re-considered by means of the second 
variation of entropy and free-energy. In contrast to the previous analysis using 
the old Tsallis formalism, the second variation of free-energy shows a distinct 
thermodynamic structure. Nevertheless, the zero-eigenvalue problem in each 
case exactly recovers the marginal stability condition inferred from the specific 
heat. Thus, within the new formalism using normalized q- values, all the anal- 
yses are consistent and the stellar polytrope can be regarded as a plausible 
non-extensive thermal state. Finally, Section 5 is devoted to the conclusion 
and the discussion. 



2 Maximum entropy principle revisited 

Throughout the paper, we pursue to investigate the equilibrium property of the 
stellar self-gravitating system consisting of N particles confined in a spherical 
cavity of radius r e . For simplicity, each particle has the same mass m and 
interacts via Newton gravity only. Then, the total mass becomes M = NrriQ. 
In this situation, owing to the maximum entropy principle, we seek the most 
probable quasi-equilibrium distribution in an adiabatic treatment. That is, 
we consider the quasi-equilibrium structure as the extremum-entropy state in 
which the particles elastically bounce from the wall, keeping the mass M and 
energy E constant. 

Following the papers I and II, we treat this issue employing the mean-field ap- 
proach that the correlation between particles is smeared out and the system 
can be simply described by the one-particle distribution function f(x, v), de- 
fined in six-dimensional phase-space (x, v). In this treatment, the one-particle 
distribution is regarded as a fundamental quantity characterizing the stellar 
system. Let us denote the phase-space element as h 3 = (lov ) 3 with unit length 
l and unit velocity v and define the integral measure d e T = d 3 xd 3 v/h 3 . Re- 
garding the function f(x, v) as a fundamental statistical quantity, the energy 
and the mass are respectively expressed as follows: 



E = K + U 



/{ 



(1) 




(2) 



where the quantity $(cc) is the gravitational potential given by 
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^ X ) = -G I d e T / (3) 

•' \x — x'\ 



On the other hand, in the new framework of Tsallis' non-extensive thermo- 
statistics, all the macroscopic observables of the quasi-equilibrium system can 
be characterized by the escort distribution, but the escort distribution itself 
is not thought to be fundamental. Rather, there exists a more fundamental 
probability function p(x, v) that quantifies the phase-space structure. With a 
help of this function, the escort distribution is defined and the macroscopic ob- 
servables are expressed as the normalized g-expectation value as follows (e.g., 
Refs.[16,17]): 



escort distribution : P q (x, v) = {p( x i v )} ^ ^ 

jd 6 r{p(x,v)r 

normalized g-value : {Oi) q = J d 6 r 0{ P q (x,v) . (5) 

And based on the fundamental probability p(x, v), the Tsallis entropy is given 
by 



Jd 6 T[{p(x,v)Y-p(x,v)}. (6) 



Note that the probability p(x, v) satisfies the normalization condition: 

J d 6 T P (x,v) = 1. (7) 

To apply the above Tsallis formalism to the present problem without changing 
the definition of energy and mass (1) and (2), we identify the one-particle dis- 
tribution with the escort distribution P q , not the probability function p(x, v): 



f(x, v) — M M^l>! ; N q = J d e r {p(x, v)r (8) 



so as to satisfy the mass conservation (2). With this identification (8), the 
normalized q- values are naturally incorporated into our mean- field treatment. 
As a consequence, later analyses become almost parallel to the previous study, 
in a form-invariant manner. 

Now, adopting the relation (8), let us seek the extremum-entropy state under 
the constraints (1) and (7). As has been discussed recently, there are two ap- 
proaches that extremize the entropy under certain constraints, i.e., the method 
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developed by Tsallis, Mendes & Plastino (TMP)[16] and the optimal Lagrange 
multiplier(OLM) method by Martinez et al. [17]. Here, we specifically apply 
the TMP procedure and find the extremum entropy state. The alternative 
derivation using OLM method is presented in Appendix A. The only differ- 
ences in the final expressions between the TMP and the OLM method are the 
dependence of the quantity N q , which can be summarized in a unified form 
(see Eq.(15) with (16)). 

The variational problem in the TMP method is given by the following equa- 
tion: 



S q -a{J d *T P -l}-p{J d «T& + U)f-E} 



= 0, 



(9) 



where the variables a and (5 denote the Lagrange multipliers. The variation 
with respect to the probability p(x, v) leads to 



J <fr — (qp^ 1 -l)6p-a6p-P (^v 2 + 5f 



0. (10) 



Here, we used the fact that / d 6 T5Qf = J d 6 T$5f. The above equation in- 
cludes the variation 5f, which can be expressed with a help of the relation 
(8): 



Then, substituting the above expression into (10) becomes 



/ 



- a - 13 M q 



' -v 2 + $ - e 



N q V2 



5p = 0. (12) 



with the quantity e being 



(13) 



In arriving at the equation (12), we have exchanged the role of the variables 
(a;, v) <-> (x', v'). Since the above equation must be satisfied independently of 
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the choice of the variation dp, we obtain 

1 , .. //' / I 



—^k 1 - 1 )-"-^'!- il v2+ ®- £ ) =o ' (i4) 



which reduces to the same power-law distribution of e as has been derived in 
papers I and II. That is, even in the new formalism, the extremum-entropy 
state remains unchanged and can be described by the so-called stellar poly- 
trope. Together with the result by OLM method (Appendix A), it can be 
summarized as follows: 



$0 - -v 2 - $(ar) 



5/(1-5) 



(15) 



where we define the constants A and 



A -w, \ a[l - q) + l Jt>) • JmT^q) (16) 

Note that the quantity N' q means N' q = N q for the TMP method and N' q = 1 
for the result using OLM procedure. 

For later analysis, we define the density p{r) and the isotropic pressure P(r) 
at the radius r = \x\ as 



= 4^2* B (| ^ A [$ _ $ (r) ]V(l-5) + l/2 



and 



. /" d 3 v 1 9 „. 

p(r) = 7 W z !(x ' v) ' 



^ b (It^)^*°-*wi i/m+3/2 < 18 > 

Here, the function £>(a, 6) denotes the beta function. These two equations lead 
to the following polytropic relation: 

P{r)=K n p l+1 l n {r), (19) 
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with the polytrope index n given by 
1 



11 

n =-T- + o> (2°) 



and with the dimensional constant K, 



n ■ 



n + 1 



n-3/2 - , " ~ X l n 



fr 3 V2' 27 [ a (l-q) + l N' q J JV, 

Using these quantities, the one-particle distribution can be rewritten as fol 
lows: 

1 p h 3 



.(21) 



f(x,v) 



4^5(3/2,71-1/2) {(n + l)K n pVn}3/2 

V 2 /2 



x 1- 



(n + l)X nP V 




which agrees with the previous result (see Eq.(16) with the identification (n + 
l)K n = (n - 3/2)T in paper I). 

While the resultant form of the quasi-equilibrium distribution (22) turns out 
to be invariant irrespective of the choice of the statistical averages, we should 
be aware of the two important differences between the old and the new Tsallis 
formalism, which we shall describe below. 

First, the relation between the polytrope index n and Tsallis parameter q in 
the polytropic relation (19) differs from the one obtained previously, but is re- 
lated to it through the duality transformation, q <-> 1/q (see Eq.(14) in paper 
I or Eq.(12) in paper II). This property has been first addressed in Ref.[16] in 
more general context, together with the changes in Lagrangian multiplier f3. 
The duality relation implies that all of the thermodynamic properties in the 
new formalism can also be translated into those obtained in the old formal- 
ism. As shown in section 3 and 4, this is indeed true in the system confined 
in an adiabatic wall (micro-canonical ensemble case), however, the duality of 
thermodynamic structure cannot hold in the system surrounded by a thermal 
bath(canonical ensemble case). 

Second, notice that the quasi-equilibrium distribution (15) with (16) contains 
the new quantities N q and e, which implicitly depend on the distribution 
function itself. In marked contrast to the result in old Tsallis formalism, this 
fact gives rise to the non-trivial thermodynamic relations as follows. Using the 
definitions of density and pressure (17) and (18), the quantity e becomes 
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£ = J4 / < PxP{.x) + J d 3 x p(x) $(x) | 

= if {I J d3xP{x) ~ J d3xp{x) [$o " ® {x)] } + $o - 

Further using the relation $ — $(x) = (n + l)(P/p) from (17) and (18) and 
substituting the equation (16) into the above expression, the variable e is 
cancelled and the equation reduces to 

^ = Jd'xP(x). (23) 



As for the dependence of N q , the normalization condition (7) implies that 



N„ 



'fd 6 T {f(x,v)}^ 



Substituting the distribution function (22) into the above equation and inte- 
grating over the velocity space, after some manipulation, one obtains 



N q = c n K^-^ [ jd\r 



l+l ' n (x) } 



-(n-3/2)/(n-l/2) 



(24) 



with the constant c n given by 

4y/2 7rE(3/2,n-l/2) / (n + lp 1 1/( "" 1/2) 

° n r r- , , , s 1 (n-3/2)/(n-l/2) 1 J.3 [ " ^ °' ) 

{4^^5(3/2,71+1/2)} I h ) 



Equations (23) and (24) play a crucial role in determining the thermodynamic 
temperature of stellar polytrope in section 3, as well as the onset of gravother- 
mal instability in section 4. In particular, the equation (23) yields the radius- 
mass-temperature relation characterizing the quasi-equilibrium structure of the 
system confined in a thermal bath. 

Keeping the above remarks in mind, hereafter, we will specifically focus on 
the spherically symmetric case with the polytrope index n > 3/2 (q > 0), in 
which the quasi-equilibrium distribution is at least dynamically stable (see 
Chap. 5 of Ref.[l]). In this case, despite the above detailed differences, the 
stellar quasi-equilibrium distribution can be characterized by the so-called 
Emden solutions (e.g., [21,22]) and all the physical quantities are expressed 
in terms of the homology invariant variables (u,v), which are subsequently 
used in later analysis. In Appendix B, together with some useful integral 
formulae, we summarize the relation between the Emden solution and the 
stellar polytropic distribution. 



9 



3 Thermodynamic properties of stellar polytrope 



Having established the quasi-equilibrium distribution, we now investigate the 
thermodynamic properties of stellar polytropic system. To do this, we adopt 
the same procedure as in paper II. That is, we examine the Clausius relation 
under the quasi-static variation for the new equilibrium system in section 3.1. 
In contrast to the previous analysis, the thermodynamic temperature can be 
identified with a help of (23), consistently with the recent claim based on the 
thermodynamic zero-th law. Using this temperature, in section 3.2, presence 
or absence of thermodynamic instability is discussed evaluating the specific 
heat of stellar system. 



3. 1 On the definition of thermodynamic temperature 



As usual, the thermodynamic instability in stellar system is intimately related 
to the presence of negative specific heat [5,23]. The evaluation of specific heat 
is thus necessary for clarifying the thermodynamic property. In this regard, 
the identification of temperature in stellar system is the most essential task. 

In the new framework of Tsallis' non-extensive thermostatistics, the physically 
plausible thermodynamic temperature, T p h ys , can be defined from the zero- 
th law of thermodynamics [24,25,26]. Then, the thermodynamic temperature 
in non-extensive system differs from the usual one, i.e., the inverse of the 
Lagrange multiplier, f3. Depending on the methods extremizing the entropy, 
one has 



phys 



/T 1 ; (OLM method), 

[1 + (1 - q) S q ] /T 1 ■ (TMP method). 



Translating the above result into our notation immediately yields that the 
physical temperature in stellar polytropic system is given by T p h ys = N'//3. 
However, the relation (26) should be carefully applied to the present case, since 
the verification of thermodynamic zero-th law is very difficult in stellar equi- 
librium system with long-range interaction. Furthermore, even using the new 
formalism, the energy E still keeps non-extensive due to the self-referential 
form of the potential energy (see Eqs.(l)(3)). In order to validate the use of 
the definition (26), we therefore adopt a rather simple procedure as exam- 
ined in paper II. That is, we consider the relation between heat transfer and 
entropy change in the quasi-static treatment under keeping the total mass con- 
stant. According to the definition (26), the Clausius relation is appropriately 
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modified and is expressed as follows [24,25,26]: 

dS q = -^—{l + (l-q)S q }d'Q. (27) 

-tphys 



This, in turn, implies that the modified Clausius relation can be used as a 
consistency check of the physical temperature (26) in stellar system. 

Let us first write down the entropy of extremum state: 

S q =(n- 1) (N q - 1). (28) 



/,From equation (24) with the polytropic relation (19), the quantity N q in 
equation (28) is rewritten as 



N q = c n K^ n ^ jy>z P(x)} 



-(ra-3/2)(n-l/2) 



f 1 GM 2 /T7 + 1 7/ . 1 "(n-3/2)/(n-l/2) 

= c„ 1^-^) -J- ™- ( ^ - 2^ - l) (29) 
[ n - 5 r e V v e v e J J 

with a help of the integral formula (88). Here, the variables (u e ,v e ) denote 
the homology invariants evaluated at the boundary r = r e (see definitions 
(84) (85)). Note also the fact that the constant K n is expressed in terms of 
(u e , t> e )-variables: 



k - p ' - GM 1 1 V" rmi 

Using these expressions, the variation of the quantities N q and K n respectively 
becomes 



1\ dN q dK n ( 3\ dr e 

n = n + \n 

2) N q K n \ 2J r e 

n — 3/2 GM 2 8 { n + 1 dv e ( du e dv e \\ , . 

2 — - — , (31) 



and 



n — 5 r e N' [ v e v e \ u P v P 



dK n n — 3 dr e dv e 1 du e 
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In the last line of equation (31), we have used the relation (23). Collecting 
these results, the entropy change dS q can be expressed as the variations of 
both the homology invariants (u e , v e ) and the wall radius r e as follows: 



dS, : - in--) dN q = N q 



3 dr e du e 
2 r e u e 



dv P 



n ■ 



n- 3/2 GM 2 (3 f n + 1 dv e ( du e dv e 



n — 5 r e N' 

e q 



(33) 



In the above equation, the term GM 2 f3 / (r e N' q ) can be factorized out using the 
relation (23) and we have 



1 1 



3 1 / Ue _ n+1 | A rfr e 
2 n — 5 V f „ v P J r P 



1 



+ — -— ^4 In - -1 u e + 2v e - 2(n + 1)) ^ 
n — 5 2f P I V 2/ J tt,, 



+ {6w e + 2m; e -3(n+l)} 



(34) 



Further notice the fact that the variation of homology invariants is expressed as 
the variation of dimensionless quantity d^ e (see Eq.(86)). The entropy change 
is finally reduced to the following expression: 



^ N GM2(3 



x <^ Au e + 2u e v e - 



3 1 



2 n — 5 V v, 

n+1 

8 + 3 ) u e - 

n — n 



u e n + 1 \ dr e n — 2 1 
2— + 11 — 



r P n — 5 2v P 



n-2 \n-2j\ ie 



(35) 



The Clausius relation in a quasi-static variation relates the variation of entropy 
with the heat change. According to the first-law of thermodynamics, the heat 
change in a quasi-static variation is estimated as follows (see Eq. (26) in paper 
II): 



d'Q = dE + P e dV = d{-\ 



GM 2 



+ Anr 2 P e dr e 



GM 2 



A+ M _ £ d^d^e) 

V e J r e e d£ e £e J ' 



(36) 



where the dimensionless quantity A related to the energy and its derivative 
dX/d^e are respectively given by (see Eq.(91)): 
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A 



r e E 
GM 2 



n 



-, il 1 1 



n + 1 



+ (n 



2)^}, 



(37) 



and 



dA _ w - 2 g(u e ,v e ) _ 
e d£ e n — 5 2f e 

,(«, .) = 4»» + 2» - {8 + 3 (2±») } „ - ^ . + 3 (2±1) . (38) 



Hence, the comparison between (35) and (36) with (37) (38) leads to the mod- 
ified Clausius relation corresponding to the equation (27): 



dS q = -ji N q (dE + P e dV) = A {1 + (1 _ q ) Sq } d'Q. (39) 

q q 

Therefore, with this relation, the thermodynamic temperature can be consis- 
tently identified along the line of the argument in Ref.[24] and the plausible 
physical temperature is now 

N' 

?Vs = -f- (40) 



3.2 Thermodynamic instability arising from the negative specific heat 



Once adopting the definition (40), we immediately obtain the radius-mass- 
temperature relation in terms of the homology invariants as follows. Equation 
(23) in section 2 leads to 

/■,,,, 1 GM 2 / u e n+1 \ 

T^ = j^P( x) = -——(2^- — + l) (41) 

from (88). Thus, defining the dimensionless quantity rj = GM 2 /(r e T phys ), we 
obtain 

GjfL = (n-5)«. 
r e T phys n + 1 - 2u e - v e 



This is in marked contrast to the result using the standard linear mean (c.f. 
Eq.(31) of paper II). While the radius-mass-temperature relation in previous 
paper includes the residual dimensional parameter h = (loV ) 3 , the expression 
(42) has no such parameter dependence and is quite similar to the result in 
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Boltzmann-Gibbs case (e.g., Eq.(29) in Ref.[5] or Eq.(25) in Ref.[27])). The 
nice form of the radius-mass-temperature relation implies that the specific heat 
can be determined independently of the residual parameter h, which is indeed 
a desirable property for a rigid theoretical prediction without any uncertainty. 

By definition, the specific heat at constant volume C v is given by 

(f) 



The numerator and the denominator in the last expression are respectively 
rewritten as follows: 



'dE\ _ GM 2 dX _ GM 2 n — 2 g(u e , v e ) 



dU P r e d£ e r e n - 5 2u e £ e 



and 



f dT phys \ _GM 2 d _ GM 2 1 k(u e ,v e ) 
\ d£ ) e r e d£ e r e n-5 u e £ e 

In the above equations, the functions g(u e ,v e ) is already given by (38) and 
the function k(u e ,v e ) is expressed as 



k(u e , v e ) = Au 2 + 2u e v e — (n + 9)u e — v e + n + 1, (46) 

from (86). Thus, the expression of specific heat can be reduced to the following 
simple form: 

n-2 g(u e ,v e ) 



Provided the homology invariants at the boundary, the specific heat in the 
new formalism is uniquely determined, free from the residual parameter h 
(c.f. Eq.(42) in paper II). In this sense, the result (47) can be regarded as a 
successful outcome of the new Tsallis formalism using the normalized q- values. 

Now, we focus on the thermodynamic instability inferred from the qualitative 
behavior of specific heat. Recall that the thermodynamic instability appears 
when the specific heat of the system changes its sign. From (47), we readily 
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expect the two possibilities. One is the case when the function g(u e , v e ) changes 
its sign. In this case, the condition for marginal stability, C v = 0, becomes 

g(u e ,v e ) = 0. (48) 



The other cases appear when the sign of the function k(u e ,v e ) is changed. In 
this case, the marginal stability leads to the divergent behavior, C v — > ±00 
and the condition is given by 

k{u e ,v e ) = 0. (49) 



To see how and when the instability develops, in Fig. 1, we plot a family of 
Emden solutions in the (77, A)-plane. Since the dimensionless parameters A and 
7] are respectively proportional to — E and T~^ ys , the signature of the specific 
heat can be easily deduced from the slope of the curve. Note that each point 
along the trajectory represents an Emden solution evaluated at the different 
value of the radius r e . /,From the boundary condition (83), all the trajectories 
start from (rj, A) = (0, —00), corresponding to the limit r e — > 0. As increasing 
the radius, the trajectories first move to the upper-right direction as marked by 
the arrow, and they suddenly change their direction to upper-left. This means 
that the divergent behavior of specific heat eventually appears and beyond 
that point, the signature of specific heat changes from positive to negative. 
That is, the potential energy dominates the kinetic energy(A > 0) and the 
quasi-equilibrium state ceases to exist for a system in contact with a thermal 
bath (paper II). Notice that even in this case, the stable quasi-equilibrium 
state still exists for a system surrounded by an adiabatic wall. On the other 
hand, for more larger radius, while the curves with index n < 5 abruptly 
terminate, the trajectories with n > 5 next reach at another critical point 
drj/dX = 0, i.e., C v = 0. Further, they progressively change their direction 
and finally spiral around a fixed point. The appearance of the critical point 
Cv = is explained as follows. While the inner part of the system keeps the 
specific heat negative, the outer part seems to have positive one. Thus, the 
heat current from inner to outer part causes the raise of the temperature at 
both parts. For a system with the radius r e smaller than certain critical value, 
the amount of the heat capacity at outer part is small so that the outer part 
easily catches up with the increase of the inner-part temperature. As increasing 
r e , the fraction of the outer normal part grows up and it eventually balances 
with the inner gravothermal part. Thus, beyond the point characterized by the 
condition Cy = 0, no thermal balance is attainable and the system becomes 
gravothermally unstable. This is even true in the system surrounded by an 
adiabatic wall. 

/,From these discussions, one can immediately verify that the condition (48) 
represents the marginal stability for a system confined in an adiabatic wall, 
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which exactly coincides with the previous result using standard linear means. 
Hence, one concludes that the quasi-equilibrium structure obtained from the 
new formalism does not alter the thermal properties in micro-canonical en- 
semble case. On the other hand, the condition (48) indicates the onset of 
thermodynamic instability for a system in contact with a thermal bath (i.e., 
canonical ensemble case), which significantly differs from the results in old 
Tsallis formalism. Of course, this might be a natural consequence of the dif- 
ferent choice of statistical average, leading to the different definition T phys , 
however, the appearance of instability in present case might seem somewhat 
curious. While the previous results indicate the unstable state at the indices 
n > 3, consistent with the suggestion by Chavanis [28], Fig. 1 implies that 
thermodynamic instability appears for stellar polytrope with any value of the 
index n. One might worry about whether the present results rigorously match 
the stability analysis from the variational problem. In next section, to check 
the consistency of new Tsallis formalism, we develop the stability analysis 
based on the second variation of entropy and free-energy. 

Finally, in Fig. 2, varying the radius r e , the specific heat C v is plotted as 
a function of density contrast, D = p c /p e for typical polytrope indices with 
n > 3/2. Clearly, the critical point \C V \ — > oo marked by crosses exists in each 
case, while the marginal stability C v = only appears when n > 5 at a certain 
high density contrast (c.f. Fig. 3 in paper II). The numerical values indicated 
by arrows represent the critical values D crit evaluated at the point C v = 0, 
which are the same results as in Table 1 of paper I (see also Fig. 5). These 
behaviors can also be deduced from the energy-radius-mass relation and the 
radius-mass-temperature relation. In Figs. 3 and 4, using the expressions (42) 
and (37), the dimensionless values A and r] for various polytrope indices are 
evaluated and plotted as a function of density contrast, respectively. 3 The 
Emden trajectories in (rj, A) plane (see Fig.l) state that the marginal stability 
for the system in contact with a thermal bath, |C V | — > oo implies the first 
turning point dr)/d£ e = 0, or equivalently drj/dD = 0. Similarly, the marginal 
stability C v = represents the condition dX/dD = 0. /,From Figs. 3 and 4, we 
readily estimate the critical density contrast at the first turning points in each 
case, which exactly coincide with the points marked by the arrows and the 
crosses, respectively. Consistently, the turning point disappears when n < 5 
in Fig. 3, while it does always exist independently of the polytrope index in 
Fig. 4. 



3 Fig. 3 is essentially the same result as in Fig. 2 of paper I. 
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4 Stability /instability criteria from the variational problems 



Previous section reveals the existence of two types of thermodynamic insta- 
bility. Then, the marginal stability conditions (48) and (49) are obtained for a 
system confined in an adiabatic wall and for a system in contact with a thermal 
bath, respectively. In this section, in order to check the consistency of these 
results, we reconsider the marginal stability criteria based on the variational 
problem. 

According to the maximum entropy principle, the stable quasi-equilibrium 
state for a system confined in an adiabatic wall is only possible when the second 
variation of entropy around the extremum state is negative, i.e., 5 2 S q < 0. 
Similarly, the stable quasi-equilibrium distribution surrounded by a thermal 
bath should have minimum free-energy, indicating the positive value of the 
second variation of free-energy, i.e., S 2 F q > 0. Thus, the condition 5 2 S q = or 
5 2 F q = readily implies the marginal stability in each case. 

To begin with, let us write down the entropy of quasi-equilibrium state. ^From 
(24) and (28), we have 

«.= (»- I) {c (T^ m - m ^)- ( - 3/2)/< - 1/2) - 1 } . (50) 

Here, just for convenience, we introduced the quantity W: 

W = J d 3 x p 1+1/n . (51) 



Note that in deriving the equation (50), we used the relation T phys = K n W 
from (23). 

Using the above expression (50), we compute the variation of entropy up to the 
second order terms. To be specific, we consider the density perturbation 5p(r) 
around the quasi-equilibrium configuration, under keeping the total mass M 
and the radius r e constant. Then the variation of entropy becomes 



n 



5N a 



2, 

21 



SW 3 ra -2 {5T r 



n- 



phys 



phys W An- 1/2 

2 



phys 



n(n- 1/4) 5W\ 3 n 5T phys 5W 



n- 1/2 



W 



2n- 1/2 T t 



phys 



w 



(52) 
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As for the variation 8W, we have 



sw = s[j d 3 x P w) = !L±1 / d 3 x 1 5 P + -L M } „v» 



(53) 



Also, we write down the variation of energy: 



5E = 5 (- Jd 3 x P + i J d 3 x p$ 

= ^ST phyB + lJd 3 x (2<Mp + 5p5<S>), (54) 

where we used the fact that / d 3 x p5$ = / d 3 x &5p. 

Below, we separately analyze the variational problem in each case. 

4-1 Stability/instability criterion from the second variation of entropy 



Let us first consider the stability/instability condition for a system confined 
in an adiabatic wall. In this case, the conservation of total energy E is always 
guaranteed and we further put another constraint 8E = 0, which yields 



ST phys = ~ f d 3 x(2$8p + 5p8$) 



from (54). Substituting the above equation into (52), a straightforward calcu- 
lation leads to the variation of entropy up to the second order, summarized as 
follows: 



5S q = 5 w S q + 5 {2) S q ; 



5^S q = - A. J fa I $ + ( n + i) ^ p Vn {x) J ^ (55) 



8 {2) S q = -N q 



+ 



3T p 2 hys n - 1/2 



J \2T p 

(^J d 3 x<$5p 



phys 

1 n-2 < - ^ 2 



2W n 
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1 (n+ l) 2 (n- 1/4) 
VP- n(n- 1/2) 



1 



n + 1 



phys 



— (J<PxpV n 6p}(f<Px*6p 



(56) 



The resultant form of the second variation of entropy S^S q seems rather 
complicated, however, recalling the fact that the background solution of stellar 
polytropic distribution always satisfies the condition S^S q = 0, the last three 
terms in right-hand-side of equation (56) can be rewritten in more compact 
form. After some algebra, we obtain 



6^S q = -N q 



J \2T P 



W 



n 



„!,;,;, 2W n 

3 n + 1 T phys - - 



3Tp hys n 3/2 



j d 3 x$ + -- 



n 



W 



Sp 



(57) 



Apart from the over-all positive constant, the above expression is indeed the 
same equation as previously obtained in paper I (see Eq.(36) in paper I, with 
the identification T p ^ ys /W = {(n — 3/2)/(n + 1)}T). Thus, just following 
the same calculation as in paper I, the stability/instability criterion from 
5^S q = can be obtained, from which one can rigorously prove that the 
marginal stability condition for a system confined within an adiabatic wall 
exactly reproduces the equation (48). As a result, for certain critical values of 
dimensionless energy and density, A cr i t and -D cr it, the onset of the gravothermal 
instability appears when the stellar polytrope with index n > 5 (see Figs. 3 
and 5). 



4-2 Stability/instability criterion from the second variation of free- energy 



Now, turn to focus on the second variation of free-energy. In paper II, usual 
definition of free-energy 

F' q = E - P~ l S q (58) 



has been used to investigate the marginal stability for a system surrounded by 
a thermal bath. As several authors recently pointed out, the usual definition 
of free-energy becomes inadequate when the physical temperature cannot be 
identified with the inverse of Lagrange multiplier (3. According to the modifi- 
cation of Clausius relation (27), the free-energy must be generalized as follows 
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[24,25]: 

F q = E — T phys In {1 + (1 - q)S q } . (59) 

Using this generalized free-energy, we consider the variation up to the second 
order, keeping the thermodynamic temperature T p h ys constant, not the total 
energy E: 

1 [5N a 1 (5N g \ 2 \ 

Substituting the equations (52), (53) and (54) into the above equation and 
keeping the terms up to the second order, after some manipulation, one obtains 

5F q = 5^F q + 5^F q ■ 



5 {1) F q = J d 3 xl^{x) + (n + l)^^p^ n {x)^Sp, (60) 

8®F q = \\Sx fa* + ^±1 P ^-\8pf) 

Note that apart from the over-all factor, the first variation 5^F q exactly 
coincides with the first variation of entropy, 5 <yl ^S q , which ensures the fact 
that extremum state of generalized free-energy is the same quasi-equilibrium 
distribution as obtained from the maximum entropy principle. Using the fact 
S^F q = 0, equation (60) is now rewritten with 

-^(/^H 2 - (62) 

Obviously, the above equation differs from the old result based on the usual 
definition of free-energy (58) (c.f. Eq.(51) in paper II). This fact readily implies 
that the thermodynamic stability has nothing to do with the dynamical sta- 
bility, which apparently contradicts with the recent claim by Chavanis [27,28]. 
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Thus, in present case using normalized g-values, it seems non-trivial whether 
the above equation indeed leads to the marginal stability condition (49). 



In what follows, restricting our attention to the radial mode of density and 
potential perturbations, we discuss the existence or the absence of perturba- 
tion mode satisfying the equation 5^F q = 0. Following the papers I and II, 
we introduce the new perturbed quantity: 

1 dQ(r) . . 

6p{r) = A^ — - (63) 



Then the mass conservation SM = implies the boundary condition Q(0) = 
Q(r e ) = 0. Substituting (63) into (62) and repeating the integration by parts, 
the second variation of generalized free-energy is reduced to the following 
quadratic form: 

S (2) F g = --J dn J dr 2 Q(r 1 )K(r 1 ,r 2 )Q(r 2 ), (64) 
o o 



where K(ri,r 2 ) stands for the operator given by 



T y, ,'n + l d ( 1 P(n) d \ G \ 

K{ri ^ ] = \ — W-Un) p( ri ) drT J + r{ ( 6Din - r2) 



, 1 ^1)^2) (65) 
nT phys dn dr 2 



Thus, the problem reduces to the eigenvalue problem and the marginal stabil- 
ity condition just corresponds to the zero-eigenvalue problem: 



Jdr f K(r,r')Q(r f )=0, (66) 


which yields 



L[Q] 



n 



1 d 



1 



n dr \Anr 2 p p dr ) r 

1 Gm(r 



Q 



phys 



dr'^P-Q(r'). (67) 
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To obtain the solution of the above zero-eigenvalue problem, we follow the 
similar procedure in paper I. First recall the fact that the following equations 
are satisfied: 

; r , o / \ ~] n — 3 Gm(r) f . , . . n — 1 Gm(r) 
L47rr 3 pr = ^, Lmr = ^. 68 



Then, these relations allow us to put the ansatz of the solution, 

Q(r) = c\ Anr 3 p(r) + c 2 m(r), (69) 

and to determine the coefficients c\ and C2 by substituting the ansatz into 
(67): 

n — 3 n — 1 

ci + c 2 = A n , 70 

n n 



where we define 
A 



l - Ur^ Q(r). (71) 



phys Q 



In addition to the above condition, we must further consider the boundary 
conditions Q(0) = Q(r e ) = 0. While the condition Q(0) = is automatically 
fulfilled, the remaining condition Q(r e ) = requires 

d 4vrr 3 p e + c 2 M = 0. (72) 



Hence, from (70) and (72), the coefficients c\ and c 2 are determined and are 
expressed in terms of homology invariants: 

n A n nu e A n 

c i = o 7 T\ — ' C2 = q 1 r\ — • t 73 ) 

n — 6 — (n — 1) u e n — 6 — (n — 1) u e 



As noted in paper I, the above coefficients are not specified completely because 
of the quantity A n in the coefficients, which depends on the solution (69) itself. 
To eliminate this self-referential structure, we directly evaluate the non-local 
term A n , leading to the consistency condition for the existence of solution (69). 
Substitution of the solution (69) with the coefficients (73) into the expression 
(71) yields 



l = _ T phys / j_ Gm(r 



n — 3 — (n — l)u e 



J dr — |47rr 3 p(r) — u e m(r)| 

o 



22 



In the above equation, the integrals in the right- hand- side can be evaluated 
using the formulae (89) and (90) listed in Appendix B. Further, the radius- 
mass-temperature relation (42) eliminates the dependence of T phys , leading 
to 



2u e + v e — (n + 1) n — 3 — (n — l)u e 



x 



[n - 5K + {2u e - 1) {3 + (n + 1) (—--)} 



After some algebra, the above equation is rewritten with the following quadratic 
form: 



= 4-Ug + 2u e v e — (n + 9)u e — v e + n + 1 
= k(u e ,v e ), (74) 

which coincides with the marginally stability criterion (49) derived from the 
condition for specific heat, C v — > 00. 

Therefore, unlike the first impression, we reach at the fully satisfactory con- 
clusion that the thermodynamic instability inferred from the specific heat is 
consistently explained from the variational problems. In other words, the ap- 
plication of new Tsallis formalism to the stability of stellar quasi-equilibrium 
system reveals a consistent thermodynamic structure of self-gravitating sys- 
tem, as well as the existence of thermodynamic instability. 



5 Discussion & conclusion 



In this paper, we revisited the issues on the thermodynamic properties of stel- 
lar self-gravitating system arising from the Tsallis entropy, with a particular 
emphasis to the standard framework using the normalized (/-values. It then 
turns out that the new extremum-entropy state essentially remains unchanged 
from the previous study and is characterized by the stellar polytrope, although 
the distribution function shows several distinct properties. Taking these facts 
carefully, the thermodynamic temperature of the extremum state was identi- 
fied through the modified Clausius relation and the specific heat was evaluated 
explicitly. The detailed discussion on the behavior of specific heat finally leads 
to the conclusion that the stability of the system surrounded by a thermal 
wall (canonical case) is drastically changed from previous result, while the on- 
set of gravothermal instability remains unchanged for a system confined in an 
adiabatic wall(micro-canonical case). The existence of these thermodynamic 
instabilities can also be deduced from the variation of entropy and free-energy 



23 



rigorously. As a result, above the certain critical values of A or D, the ther- 
modynamic instability appears at n > 5 for a system confined in an adiabatic 
wall. As for the system in contact with a thermal bath, the onset of ther- 
modynamic instability appears at (^ cr i t , -D CT it) with any value of n. In Figs. 
5 and 6, for the sake of the completeness, the critical values evaluated from 
the conditions (48) and (49) are summarized, respectively. Note also the fact 
that the critical values in Fig. 5 exactly coincide with the previous results (c.f. 
Table 1 in paper I). 

The one noticeable point in the present result is that the macroscopic relation 
such as the radius-mass-temperature relation or the specific heat as a func- 
tion of density contrast can be solely determined from the Emden solutions 
without any uncertainty. Indeed, previous study using the standard linear 
means seriously suffer from the residual dimensional parameter h = (IqVq) 3 , 
which must be practically disappeared from the macroscopic description. In 
paper II, the origin of this residual dependence has been addressed in con- 
nection with the non-extensivity of the entropy. Although the scaling relation 
appeared in the old results is simply deduced from the asymptotic behavior 
of the Emden solutions, the explicit /i-dependence itself has originated from 
the radius-mass-temperature relation. By contrast, in present case, the radius- 
mass-temperature relation was derived from the non-trivial relation (23), in 
which no such /i-dependence appears. The resultant specific heat is thus ob- 
tained free from the residual dependence, which is a natural outcome of the 
new framework using the normalized g-values. Therefore, it seems likely that 
the new formalism provides a better characterization for non-extensive meta- 
equilibrium state. 

In fact, for several application of the new formalism, the equation of state 
for non-extensive system has been found to be similar to the result from the 
ordinary extensive thermodynamics. For instance, the equation of state for 
the classical gas model without interaction reduces to the ideal gas sate even 
in the power-law nature of velocity distribution [24]. Indeed, this fact can be 
clearly seen in our case neglecting the gravity. Recalling that both the density 
and the pressure become homogeneous in the limit G — > 0, the equation (41) 
immediately reduces to the relation T phys = PV, i.e., P oc pT phys , where V 
denotes the volume of the system. This result apparently seems contradiction 
with the polytropic relation (19), however, it turns out that the dimensional 
constant K n immediately yields the relation K n = Pj p 1+l / n from the definition 
(21) and thereby the polytropic relation makes no sense in the limit G — > 
0. Thus, the stellar polytropic system using normalized q- values successfully 
recovers the ideal gas limit, in contrast to the old results (see Appendix B 
in paper II). In other words, if the interaction is turned on, the equation of 
state for the ideal gas is no longer valid and instead the polytropic relation 
holds. In Sec. 2, the polytropic relation was derived without assuming any 
specific choice of the gravitational potential. This in turn suggests that at 
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least in the mean- field treatment, the polytropic relation is a common feature 
in presence of long-range interaction and it generally holds for non-extensive 
quasi-equilibrium system. 

Finally, we have firmly confirmed that the stellar polytropic system can be 
consistently characterized as a plausible meta-equilibrium state in the new 
framework of the non-extensive thermostatistics, as well as in the old for- 
malism. At present, however, one cannot rigorously discriminate the correct 
and the applicable formalism among them (there might be another possibility 
that both the formalism become indeed correct depending on the situations). 
While the new framework by means of the normalized g-values has theoret- 
ically desirable properties as mentioned above, the results in old formalism 
might provide an interesting connection with dynamical instability of gaseous 
system, as has been suggested by Chavanis [27,28]. In paper II, the thermody- 
namic instability of stellar polytropic system in the canonical ensemble case 
is shown to exist at the polytrope indices n > 3, which is exactly the same 
condition as derived from the gaseous system. On the other hand, in the new 
formalism, there is no such correspondence, since the instability is completely 
different from that of the gaseous system (see Sec. 4. 2). This means that the 
choice of the statistical average is crucial for the thermodynamic properties 
and the physical reason for this discrepancy should be further clarified. In 
any case, in order to pursue the physical reality of the non-extensive thermo- 
statistics, the limitation of thermodynamical approach is now apparent and 
the detailed kinematical study based on the Boltzmann or the Fokker-Planck 
equation provides an deep insight into the thermal transport property. In the 
light of this, the long-term stellar dynamical evolution by iV-body simulation 
also becomes useful and the analysis is now in progress. The results will be 
presented elsewhere. 

We are grateful to Prof. S. Abe for invaluable discussions and comments. This 
work is supported in part by the Grant-in-Aid for Scientific Research of Japan 
Society for the Promotion of Science (No. 1470157). 
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Appendix A: Extremum-state of Tsallis entropy from the OLM 
method 



Here, owing to the maximum Tsallis entropy principle, we derive an ex- 
tremum state of the stellar quasi-equilibrium distribution by means of the 
OLM method. According to Ref.[17], with a help of the relation (8), the vari- 
ational problem that extremizes the entropy under the energy constraint be- 
comes 



S q - a {Jd 6 rp- lj -0 J d 6 r^M (^v 2 + ^ - £ j p q 



0.(75) 



Then the variation with respect to probability function p(x, v) leads to 



/ 



Sp 



-\ps(Ng |d 6 T$/)=0. (76) 
In the above equation, the last term in the left- hand- side is rewritten with 

5 (N q J d 6 r $ /) = 5N q J d 6 r $ / + 2 N q J d e r $5/ 

= qj rf 6 r|2M$- J d 6 r' $(x')f(x', «')} P"' 1 Sp, (77) 

with a help of the equation (11). In the last line, we have exchanged the role 
of the variables (x,v) <-> (x',v'). Substituting (77) into (76), we arrive at 



J d 6 T ^(qp^-l) -a-(3Mqp q - l (^v 2 + §-e^ 



5p = 0, 



where the quantity e is given by (13). Recalling the fact that the above equa- 
tion holds independently of the choice of the variation 5p, we obtain 



- — y (qp^ 1 - l) - a- (3Mqp q ~ l {^v 2 + $ - = 0, 

which leads to the power-law distribution. The resultant expressions for the 
one-particle distribution function (escort distribution) is summarized in equa- 
tions (15) and (16), with the identification N' q = 1. 
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Appendix B: Stellar polytropic system characterized by Emden so- 
lutions 



In this appendix, we briefly describe the equilibrium configuration of stellar 
polytrope in the spherically symmetric case, together with some useful formu- 
lae which has been used in the main analysis of section 3 and 4. 

First notice that the one-particle distribution function (15) does not yet com- 
pletely specify the equilibrium configuration, due to the presence of gravita- 
tional potential which implicitly depends on the distribution function itself. 
Hence, we need to specify the gravitational potential or density profile. ^From 
the gravitational potential (3), it reads 



Combining the above equation with (17), we obtain the ordinary differential 
equation for $. Alternatively, a set of equations which represent the hydro- 
static equilibrium are derived using (78), (17) and (18): 



= -^p P (r), (79) 
dr r z 

dm(r) 



dr 



= Anp(r)r 2 . (80) 



The quantity m{r) denotes the mass evaluated at the radius r inside the wall. 
Denoting the central density and pressure by p c and P c , we then introduce the 
dimensionless quantities: 



(n + l)P c 



1/2 



which yields the following ordinary differential equation: 

6" + -& + 9 n = 0, (82) 



where prime denotes the derivative with respect to £. The quantities p c and P c 
in (81) are the density and the pressure at r = 0, respectively. To obtain the 
physically relevant solution of (82), we put the following boundary condition: 

9(0) = 1, 0'(O) = 0. (83) 
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A family of solutions satisfying (83) is referred to as the Emden solution, which 
is well-known in the subject of stellar structure (e.g., see Chap. IV of Ref.[21]). 

To characterize the equilibrium properties of Emden solutions, it is convenient 
to introduce the following set of variables, referred to as homology invariants 
[21,22]: 



d\nm(r) Anr 3 p(r) £9 r 
d\ar m{r) 9' 



dlnP(r) _ p(r) ft*) 



dlnr P(r) r 6 

which reduce the degree of equation (82) from two to one. The derivative of 
these variables with respect to £ becomes 

du ( n \ u dv ( 1 \ v 

_= 3-u -v -, -t-= [-1 + U+ -v -. 86 

d£ V n + 1 / £ rf£V n + 1 / £ 



Equations (82) can thus be re-written with 

« (it) (n + l)(u — 1) + v 
v du (n+ 1)(3 — u) — nv 



(87) 



The corresponding boundary condition to (83) becomes (u,v) = (3,0). 

Now, utilizing the above homology invariants, we list some useful formulae 
which can be derived from the hydrostatic equations (79) and (80) (see also 
Appendix A in paper I): 



IdrAnr 2 P(r) = — Lnr 3 ^ - (n + 1)^ + — 1 

J n-5 [ p e r e J 



n — 5 r, 



CM 2 / 2 « e Ti+l | 1( 

V P V P 



(88) 



J dr Gm ^ 47rr 3 p(r) = -Anr 3 P e + 3 ^ dr Anr 2 P(r), 



n + 1 CM 2 



71 — 5 7V 



3 3 
— + 



71+1 



(89) 
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dr 



Gm 2 (r) GM 2 



l e 

8nr 3 e P e + 6j dr 4nr 2 P(r), 



n+1 GM 2 



n 



o 



(90) 



Finally, using these formulae, we evaluate the total energy of the stellar system 
in terms of the variables at the edge r e : 



/. : K I -- I dr i7vr 2 P<r) 
o 



J drAnr 2 P(r) - J 



d,t 



Gm(r) dm 
r dr 



1 



n — 5 



3 f GM 2 . „ . MP e } . ^ A 3 „ 

(n+1) e - \ + (n-2)Airr 3 e P ( 

2 [ r e p e I 



which can be re-expressed in terms of the homology invariants: 



E = 



1 GM 2 



n — 5 r e 



3 f n + 1 

2 I v fi 



(91) 
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Fig. 1. Trajectory of Emden solutions in (77, A)-plane. 
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Fig. 2. Specific heat as a function of density contrast D{= p c /p e ) for various poly- 
trope indices. 
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Fig. 3. Energy-radius-mass relation as a function of density contrast. 
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Fig. 4. Radius-mass-temperature relation as a function of density contrast. 
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